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Abstract. We introduce a new Macaulay 2 package, SimplicialDecomposability, which 
works in conjunction with the extant package SimplicialComplexes in order to compute a 
shelling order, if one exists, of a specified simplicial complex. Further, methods for de- 
termining vertex-decomposability are implemented, along with methods for determining 
■ fc-decomposability. 

hj ; 

Introduction. A simplicial complex A on a finite vertex set V is a set of subsets of V closed 
under inclusion. Elements a G A are called faces and the dimension of a is #cr — 1. The 
dimension of A is max dim ex. The / -vector of a A, where d = dim A + 1, is the (d+ l)-tuple 
. . ■ , /<2-i), where /, is the number of faces of dimension i in A. Using this, the h-vector 
of A is the d + 1-tuple (h , ...,h d ) given by h, = £f =0 (- 1)^01*) £-1 for < j < d. 

The Stanley- Reisner ideal is the ideal /(A) generated by the minimal non- faces of A 
and the Stanley- Reisner ring is the ring K[A] = K[V]/I(A), for a given field K. Thus 
the Stanley-Reisner ideals of complexes on a given vertex set V are exactly the squarefree 
monomial ideals in lf[V]. Using relations between the complex and the ideal, one can use 
tools from both algebra and combinatorics to study properties of both. For example, the 
/i-vector of a complex A is the coefficient-vector of the numerator of the Hilbert series of 
K[A}. 

The package SimplicialComplexes by Sorin Popescu, Gregory G. Smith, and Mike Stillman 
already implements many methods for simplicial complexes in Macaulay 2 |M2j . a software 
system designed to aid in research of commutative algebra and algebraic geometry. We intro- 
duce a new package, SimplicialDecomposability, for Macaulay 2 which provides several new 
methods for testing various forms of decomposability for simplicial complexes. Particularly, 
the package implements methods for testing shellability and vertex-decomposability. 

Shellability. If a simplicial complex has one facet, say a, then it is a simplex and is denoted 
2° ' . Let A be a simplicial complex which has equi-dimensional facets, i.e., is pure. Then by 
tj Definition III.2.1], A is shellable if its facets can be ordered o~x, . . . ,a n such that 

i i—l 



has a unique minimal element for 2 < i < n, such an ordering is called a shelling order. 

See |BW-14 Definition 2.1] for the definition of non-pure shellability, which is implemented 
in the package for non-pure complexes. 
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Shellability is of interest because it implies a number of nice properties. In particular, if 
a pure simplicial complex is shellable, then its Stanley-Reisner ring is Cohen- Macaulay over 
every field (StJ Theorem III. 2. 5]. Hence, its h- vector is non-negative and can be read off 
from any shelling order [StJ Theorem III. 2. 3]. Further still, the h- vectors of pure shellable 
complexes are numerically classified [StJ Theorems II. 2. 2 and II. 3. 3]. 

We recall that the Alexander dual of a simplicial complex A on vertex set V is the simplicial 
complex A v := {V\F | F ^ A}. Further, we say an ideal I = (fx, . . . , /„) has linear quotients 
if for 1 < i < n, the quotient ideal (fx, . . . , fi-\) : (fi) is generated by linear forms. 

In the following example we demonstrate |1 11 IX, Theorem 1.4(c)] which shows that a pure 
simplicial complex is shellable if and only if the Stanley Reisner ideal of the Alexander dual 
has linear quotients. We begin by constructing the polynomial ring R — Q[a, b, c, d, e, /, g] 
and a simplicial complex D, which we verify is pure. Note that loading the package 
SimplicialDecomposability automatically loads the package SimplicialComplexes. 

11 : needsPackage "SimplicialDecomposability"; 

12 : R = QQ[a. .g] ; 

13 : D = simplicialComplex monomial Ideal {a*b,a*c,b*c,c*d,d*e,d*f ,f *g}; 

14 : isPure D 
o4 = true 

We can recover the sequence of linear quotients directly from a shelling order. We recall 
that a pure simplicial complex A is shellable if there is an order of the facets Fx , . . . , F n such 
that for < j < i there exists an x G Fi\Fj and a < k < i such that Fi\Fk = {x}. The set 
of vertices associated to each % in the preceding statement generate the linear quotient order 
of /(A v ) with respect to the given shelling order (see the proof of [HHZl Theorem 1.4(c)]). 

15 : — find the linear quotients from a shelling order 

linearQuotients = -> ( 
for i from 1 to #0-1 list 

unique flatten for j from to i-1 list ( 
ImJ = set support 0_i - set support 0_j ; 
for k from to i - 1 list ( 

ImK = set support 0_i - set support 0_k; 
if #ImK == 1 and isSubset(ImK, ImJ) then 
first toList ImK else continue))); 

We generate a shelling order Ox of D with the method shellingOrder. This method 
attempts to build up a shelling order of D recursively using a depth-first search, adding one 
facet at a time. We note that in the non-pure case, the method only searches the remaining 
facets of largest dimension. 

16 : 01 = shellingOrder D 

06 = {c*e*g, b*e*g, a*e*g, b*d*g, a*d*g, c*e*f , b*e*f, a*e*f} 

06 : List 

17 : linearQuotients 01 

07 = {{b}, {a}, {d}, {d, a}, {f}, {f, b}, {f, a}} 
o7 : List 

It is sometimes beneficial to have more than one shelling order for a given simplicial 
complex. We can use the option Random with the method shellingOrder to first apply a 
random permutation to the facets before preceding with the recursion. 
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18 : 02 = shellingOrder (D, Random => true) 

08 = -[b*e*g, a*e*g, a*d*g, b*e*f, c*e*f , a*e*f , c*e*g, b*d*g} 

08 : List 

19 : linearQuotients 02 

09 = {{a}, {d}, {f}, {c}, {f, a}, {c, g}, {d, b}} 

09 : List 

Alternately, we may use the option Permutation with the method shellingOrder to force 
a given permutation on the facets before preceding with the recursion. 

110 : 03 = shellingOrder (D, Permutation => {3,2,1,0,4,5,6,7}) 

010 = {b*d*g, b*e*g, a*e*g, c*e*g, a*d*g, c*e*f , b*e*f, a*e*f} 

010 : List 

111 : linearQuotients 03 

011 = {{e}, {a}, {c}, {a, d}, {f}, {f, b}, {f, a}} 

011 : List 

Thus we now have multiple distinct linear quotient orders associated to the ideal I(D V ), 
each coming from a distinct shelling order of D. 

Vertex-decomposability. Let A be a pure simplicial complex and o a face of A. Then 
the link and face deletion of A by a are the simplicial complexes 

link A a := {t E A \ a f) r = ®, a U r E A} and del A cr := {r G A | a <£ r}. 

Then |PB[ Definition 2.1] defines A to be vertex-decomposable if either A is a simplex or 
there exists a vertex x G A, called a shedding vertex, such that link A x and del A x are 
vertex-decomposable. 

See |BW-2t Definition 11.1] for the definition of non-pure vertex-decomposability, which is 
implemented in the package for non-pure complexes. Also, see |Wot Definitions 3.1 and 3.6] 
for the generalisation of vertex-decomposability, called fc-decomposability. It is implemented 
in the package with the methods iskDecomposable and isSheddingFace. 

Being vertex-decomposable is a strong property which implies many things. A pure vertex- 
decomposable simplicial complex is shellable |PB1 Theorem 2.8] and hence has non-negative 
/i-vector jStJ Theorem III. 2. 3] and its Stanley-Reisner ring is Cohen-Macaulay JStJ Theo- 
rem III. 2. 5]. Furthermore, the h- vectors are numerically classified for vertex-decomposable 
simplicial complexes |Leet Theorem 3.5]. Moreover, the Stanley-Reisner ring of a pure 
vertex-decomposable complex is squarefree glicci, that is, in the Gorenstein liaison class 
of a complete intersection such that the even links are squarefree monomials |NR| Theo- 
rem 3.3]. 

In the following example we demonstrate that the simplicial complex D from the previous 
example is indeed squarefree glicci. We use |NRl Remark 2.4] to find a basic double link of 
1(D) to /(linkof), both in R, for some shedding vertex v of D. 

First, we verify that D is vertex-decomposable. We then find all of its shedding vertices. 

112 : isVertexDecomposable D 

012 = true 

113 : select (allFaces(D, 0), v -> isSheddingVertex(D, v)) 

013 = {a, b, c, d, f} 
ol3 : List 
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We choose the shedding vertex f of D and generate E = link/) /. Moreover, we in turn 
find its shedding vertices. 

il4 : E = link(D, f ) ; ideal E 

ol5 = ideal (a*b, a*c, b*c, d, f, g) 

015 : Ideal of R 

116 : select (allFaces (E, 0), v -> isSheddingVertex(E, v)) 

016 = {a, b, c} 
0I6 : List 

We now choose the shedding vertex c of D and generate F = links c. Notice then that 
the Stanley Reisner ideal of F is a complete intersection. 

117 : F = link(E, c) ; ideal F 
0I8 = ideal (a, b, c, d, f, g) 
0I8 : Ideal of R 

Hence, we now have the following sequence of basic double links in R which has squarefree 
terms on the even steps (the odd steps are omitted): 

Q[D] = (ab, ac, be, cd, de, df, fg) ~ Q[E] = (ab, ac, be, d, f, g) ~ Q[F] = (a, b, c, d, f, g). 
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